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ABSTRACT 

We report the spectroscopic confirmation of a sub-mm galaxy (SMG) at z = 4.547 with an estimated 
L/_r = 0.5— 2. Ox 1O 13 L . The spectra, mid-IR, and X-ray properties indicate the bolometric luminosity 
is dominated by star formation at a rate of > IOOOMq^t- -1 . Multiple, spatially separated components 
are visible in the Ly- Alpha line with an observed velocity difference of up to 380 km/sec and the 
object morphology indicates a merger. The best fit spectral energy distribution and spectral line 
indicators suggest the object is 2-8 Myr old and contains > I0 10 M Q of stellar mass. This object is a 
likely progenitor for the massive early type systems seen at z ~ 2. 

Subject headings: galaxies: evolution, galaxies: formation, galaxies: high-redshift, galaxies: interac- 
tions, galaxies: starburst, submillimeter 



1. INTRODUCTION 

The study of galaxies detected at millimeter (mm) and 
sub-mm wavelengths is one of the most rapidly develop- 
ing fields in observational astronomy. It is now known 
that a large fraction of the star formation activity is en- 
shrouded in dust, with the star formation rate (SFR) 
being directly proportional to the far infrared (FIR) lu- 
minosity of galaxies, modulo possib le contributions from 
an Active Galactic Nuclei ( AGN) (|Hughes et all [1991 . 
Surveys performed at mm wavelengths directly probe the 
FIR luminosity, and hence the amount of star formation. 
Furthermore, the shape of the galaxy spectral energy dis- 
tributions (SEDs) at rest-frame mm wavelengths results 
in a negative K-correction between 0.5 < z < 10. There- 
fore a flux limited survey is equivalent to a SFR limited 
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survey at these redshifts (jBlain et al.ll2002f ) . 

The current redshift distribution of m m galaxies peaks 
at z ~ 2, with very few galaxies at z > 3 ([Chapman et aT] 
120051: iPope et alj|2005t lAretxaga et al.ll2007f) . However, 
the small bandwidth of current mm wave spectrographs 
makes it very difficult to measure redshifts directly, and 
the low angular resolution of mm single dish imaging 
leads to multiple optical candidates for the same source. 
As a result, mm surveys have relied on high resolution 
radio data to identify the optical counterparts for subse- 
quent spectroscopic follow-up. This leaves 35 — 70% of 
the population of mm-selected galaxies at milli-Jansky 
flux levels unident i fied and potentially at higher redshift 
(|Wang et al.ll2007t lYounger et al.ll2007ft . 

In this paper we report the discovery of a mm galaxy 
with a spectroscopic redshift of z = 4.547 that appears to 
be dominated by star formation. This is the highest red- 
shift galaxy detected at mm wavelengths not associated 
with an optically bright quasar. The object was inde- 
pend ently selected as a Lyman-break galaxy (|Lee et al.l 
l2008h for s pectroscopic follow up, a mm source, and a ra- 
dio source (|Carilh' et al.ll200ct ). The source reported here 
is unusual for mm sources because it has several nearby 
optically bright counterparts which were selected as V 
band dropouts and is unusually luminous which allows 
for a radio detection. However, the confirmation of this 
object suggests the population of sources with similar ra- 
dio to mm flux ratios and optical colors may also be at 
high redshifts. 

We assume a H Q = 70, fl v — 0.7, Q m — 0.3 c osmology 
and a star formation rate integrated across a iSalpeterl 
(|1955f) IMF from 0.1-100A/© throughout this paper. 

2. DATA 

Observations at A = 1.1mm with an avera ge RMS noise 
of 1.3 mJy were obtained with the AzTEC ( Wilson et alJ 
120081) cam era at the the Jam es Clerk Maxwell Telescope 
(JCMT) (|Scott et all [2008). Additional observations 
were obtained by the MAMBO camera on the IRAM 30m 
telescope with an RMS of 0.67 mJy and a positional ac- 
curacy of < 5" (Schinnerer et al. in prep). Ground 
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Fig. 1. — Images in Bj, Ly-a, i + , ACS F814W, z + , K s , and 3.6/xm are shown for 6" X 6" (39.4kpc) around the source. The radio position 
is marked with a cross and a foreground object is marked with an arrow. 1.4GHz Radio contours are overlaid on the ACS band image in 
the right most panel. The Ly-a image is generated by subtracting the r + continuum image from the IA679 image which is centered on 
the Ly— a line. The 3.6/^m (rest frame optical) flux is centered on the radio position and diffuse i + and z + (rest frame UV) flux. A UV 
bright knot is visible down and to the left (0.6" to the SE) of the radio position and contains (73%) of the UV flux. The Ly-a emission is 
more extended than the UV emission seen in the i + band and originates from both the UV bright knot and the radio position but not the 
foreground objects. Note the UV bright knot is absent in the K s band image, indicating a very young age. 
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Fig. 2. — The 2-D spectra of the object around the Ly-a line is 
shown on the right along with the ground based z + band image 
on the left with the slit position marked in blue, the radio position 
marked with a magenta cross, and the foreground object indicated 
by a magenta arrow. The red lines are spaced at 0.4" intervals 
which is comparable to the spacial resolution of the spectra and 
the black bar in the bottom right has a length of lA in the rest 
frame. Ly-a emission is coming from all portions of the source 
and there is a clear velocity gradient along the slit. Note the Ly-a 
absorption feature at the physical position of the bright continuum 
emission and the multiple peaks in the Ly-a line. 



based optical and near infrared imaging in 22 bands, 
Hubble Space Telescope, Spitzer, and Very Large Array 
images were obtained as part of t he COSMOS survey a s 
described in ICapak et alJ (12007T). IScoville et all (2 007). 
ISanders et all (|2007f) . and ISchinnerer et alj (|2007f ) re- 
spectively. Additional deep J and K s data were obtained 
with the UKIRT, and CFHT telescopes as part of the 
COSMOS survey (McCraken et al. in prep). X-ray data 
were obtained with the Chandra X-ray observatory as 
part of the C-COSMOS program (Elvis et al. in prep). 
The respective fluxes are tabulated in Table [TJ 

The optical spectra were taken on the Keck-II telescope 
using the DEIMOS instrument with a 1" slit width, the 
830 line/mm grating (3.3A FWHM resolution), and the 
OG550 blocking filter to optimize red throughput. The 
data were collected in eight f800s exposures (4h total 
integration) under photometric conditions with 0.4-0.6" 
seeing. The object was dithered along the slit by ±3" 
between exposures to improve background subtraction. 
The data were reduced wi th a modified versio n of the 
DEEP2 DEIMOS pipeline (jMarinoni et alJl2001f >. In ad- 
dition to the standard processing this modified pipeline 
constructs and subtracts a median background and ac- 
counts for the shifts when combining the spectra. These 
additional steps remove the "ghosting" inherent to the 
8301/mm grating. 

An image of the continuum and Ly-a emission is shown 



Wavelength 


Flux (AiJy)+ 


0.2-8Kev 


< 3 X 10~ 16 ergs s~ 


150nm (FUV) 


< 0.2 


250nm (NUV) 


< 0.09 


380nm (u*) 


< 0.01 


427nm 


< 0.02 


446nm (Bj) 


< 0.01 


464nm 


0.04 ± 0.04 


478nm (g+) 


0.03 ± 0.02 


484nm 


0.04 ±0.04 


505nm 


0.09 ±0.05 


527nm 


0.08 ± 0.04 


548nm (Vj) 


0.15 ±0.02 


574nm 


0.18 ±0.05 


624nm 


0.25 ±0.06 


630nm (r+) 


0.50 ±0.03 


679nm 


1.24 ±0.06 


709nm 


1.30 ±0.07 


711nm 


1.32 ±0.13 


738nm 


1.59 ±0.09 


764nm 


1.60 ±0.05 


767nm 


1.59 ±0.09 


815nm 


1.60 ±0.09 


827nm 


1.48 ±0.09 


904nm (z+) 


1.82 ±0.10 


1.25/im (J) 


2.4 ±0.8 


2.15/im (K 3 ) 


3.7±0.5 


3.6/^m 


7.9 ±0.2 


4.5^tm 


5.8 ±0.4 


5.8^tm 


3.4 ±1.3 


S.Oum 


10 ±3.6 


24/^m 


26 ± 13* 


1.1mm 


4800 ± 1500 


1.25mm 


3410 ± 670 


20cm 


45 ± 9 



All limits arc lcr limits. A nearby bright source was modeled and 
subtracted to make this measurement, the formal upper limit is 70^A.Ty 
if the nearby source is not subtracted. 

in Figure [1] and indicates emission from all components 
of the source. The 2-D and 1-D spectra are shown in 
Figures [2] and [3] respectively. Ly— a emission is detected 
from both the compact and extended portions of the ob- 
ject with a velocity gradient across the slit. Ly-a emis- 
sion from the diffuse region is redshifted with respect to 
the compact source, and a deep Ly— a absorption feature 
is present at a redshift corresponding to other absorption 
features seen in the spectra. The dispersion of distinct 
peaks in the Ly-a emission is 380km s _1 and the line 
asymmetry indicates outflow winds of up to ~ 1800km 
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Fig. 3. — The optical spectrum of the source (b lack) is shown 
along with the composite LBG spectrum (red) from Shaplcy ct al. 
(2003). The subtracted sky spectra is shown in green for compar- 
ison and the atmospheric A and B absorption bands are marked 
with boxes. The inset panel on the right shows the region around 
the Ly-a line. Lines from the interstellar medium and stellar pho- 
tospheres such as O 1,0 II,C II,C IV,Si II,Fe II, He II, are clearly 
visible in the spectra. The presence of the Si IV 1297A and C IV 
1549A P-cgyni lines, along with the He II 1640A emission lines 
indicates the presences of both W olf-Rayet and O sta rs, placing 
the age of the burst at a few Myr l lPettini et alj|200Cf ). Note the 
absence of a common forbidden He [II] absorption feature at 1264 A 
which is generated in shocked gas, suggesting a we are not looking 
through a significant column of gas. 

s _1 , typical of a merger with heavy star formation. Sev- 
eral interstellar absorption features are clearly seen in the 
continuum yielding a best fit redshift of 4.547±0.002, but 
the dispersion in the Ly-a line suggests some components 
could be ±0.02 from the continuum redshift. 

A foreground object with an i + band flux of 0.58^iJy 
is visible to the west of our source in the B band im- 
age (see Figure [lj and its spectrum shows [O II] emis- 
sion at z — 1.41. The mass of this system is too low 
to significantly gravitationally lens our system and it is 
unlikely that the mm flux originates from this source 
because the object is outside the radio and Spitzer posi- 
tional error and has an SED which indicates little obscu- 
ration. Furthermore, if the source of the mm emission 
is at z — 1.41 the 24/im, mm, and radio flux place the 
far infrared luminosity at ~ 10 14 L© and the dust tem- 
pera ture at less than 20K, which is physically unreason- 
able (iSolomon fc Vanden Boutll2005l ; lDale fc HeToil2002t 
ICarilh fc Yunlll999D . 

3. PAN-CHROMATIC PROPERTIES AND MORPHOLOGY 

The rest-frame ultraviolet (UV) properties indicate a 
merger, typical of star-forming mm sour c es seen at z ~ 2 
(Cha pman et all 120051: iPope et all 120051 : lAretxaga et~aLl 
2007). At least two distinct components are visible in the 
ACS image, and a region of extended emission is visible 
in the ground based i + and z + images which are more 
sensitive to extended emission than the HST data. The 
rest frame UV is centered at 10:00:54.516, +2:34:35.17 
with the radio and rest frame optical emission centered 
at 10:00:54.48, +2:34:35.9. 

The majority (73%) of the rest frame UV flux orig- 
inates from a compact (< 2.7kpc) knot of emission at 
the south east corner of the object, with the remain- 
ing UV emission extended over 17.7 kpc (2.7" ), while 
the rest-frame optical (Spitzer-IRAC) flux of the source 
is centered on the diffuse UV and radio emission. The 
IA679 filter corresponds to the rest frame Ly-a line and 
shows strong emission from all components of the source 
and an extended Ly-a halo around the source (see Figure 
[]}. In addition, a significant excess of flux is measured 
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Fig. 4. — The pan-chromatic SED of the source is shown along 
with the best fit SED m odel. The UV br ight portion of the source 
is fit by a 2.5Myr old Maraston (2005) stellar population with 
1.9 X 10 s Mq of stellar mass and no extinction, while the UV 
diffuse portion is fit by a 6.5Myr old population with 9.7 X 10 9 Mg 
of stellar mass and 1.4 magnitudes of visual extinction. The 8um, 
24/jm, mm, and radio portion of the SED is fit by a 2 X 10 1 Lp, 
ICharv fc Elba3 ll200lh model (blue) or a 45K IDale'fe Heloul pOOa) 
model (red). 



TABLE 2 
Best Fit SEDs 



Emission 


Age 


Av 


Mass 


UV SFR 


x 2 


Model 


Source 


(Myr) 




(A/ ) 






Diffuse 


6.5 


1.4 


1 x 10 10 


3200 ± 550 


24.3 


M05 




7.6 


1.5 


2 x 10 10 


4200 ± 730 


24.7 


BC03 


Compact 


2.5 


0.0 


2 x 10 8 


250 ± 20 


54.5 


M05 




2.9 


0.0 


6 x 10 s 


250 ± 20 


55.7 


BC03 



in the 3.6/xm band which is centered on the Ha line at 
z = 4.547. The Ly-a image, the 3.6^m excess, and the 
optical-IR colors of the diffuse region are consistent with 
all components residing at z = 4.547. 

The possible presence of Ha in the 3.6/xm band and [O 
II] in the K s band combined with the low S/N of the J 
band and multiple components with different SEDs make 
it difficult to constrain the age and mass of this object 
with stellar models because the 4000A break strength is 
degenerate with the line ratio in some cases. To reduce 
these degeneracies Ha and [O II] line flux is added to the 
stellar models i n proportion to the unobscured UV star 
formation rate (jKennicuttl fl998) . The diffuse portion of 
the source fades rapidly between the K s and z + bands, 
indicating significant obscuration. The UV compact re- 
gion is not detected redwards of the z + band, indicating 
a very young (< 5 Myr old) stellar population and little 
dust obscuration. 

A two component fit to the total integrated light pro- 
duced poor results. The UV and optical light are spa- 
tially separated, but the best fit model places the ma- 
jority of both the UV and optical light in a single un- 
obscured > O.lGyr old population. To overcome this de- 
generacy we attempt to deconvolve the UV bright knot 
and the diffuse component of the source. Flux measure- 
ments were made in all of the ground based images for 
the UV diffuse portion of the source using a 1" diame- 
ter apertu re on the original s tacked images without PSF 
matching (|Capak et alj|2007t ). An aperture correction of 
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1.4 magnit udes, based on th e ACS m orphology, was then 
applie d. A lMarastonl (|2005f ) (M05) or lBruzual fe Charlotl 
( 2003) (BC03) single b urst model with solar metallicity, 
a Salpeter IMF and a Calzetti et al.l ()2000[ ) extinction 
law were then fit to these and the IRAC measurements. 
The model fluxes were then subtracted from the total 
flux, and a second model was fit to the remaining flux. 
The resulting fits are tabulated in Table [2] and shown in 
Figure [H 

The present data do not constrain the peak or shape 
of the FIR emission, so the total luminosity is uncertain, 
but the mm and 24/im data imply an infrared luminos- 
ity of 0.5 — 2 x 10 13 L©, and a correspondi ng star forma- 



tion rate (SFR) of 1000 - AOOOM^yr- 1 JlDale fc Heloul 
[200llCharv fc Elbazll200l iCarilli fcYunll 1999ft . The ra- 
dio flux gives a second estimate of the SFR independent 
of dust obscuration at 3700 ± IQOM^yr -1 a ssuming a 
radio spectral slope of -0.8 (|Kennicuttl Il998) . Finally, 
assuming the excess flux in the 3.6/xm band is due to 
Ha originating from the diffuse component, we derive a 
dust corrected star formation rate of 2900 ± 100M©yr _1 
(jCalzetti et al.ll2000tlKennicuttlll998ft . in good agreement 
with the UV, but lower than the mm, and radio deter- 
mined values. This last measurement is largely indepen- 
dent of the models because the Ha flux was fixed to the 
model star formation rate rather than fit as an extra pa- 
rameter. 

The bolometric luminosity of this source appears to 
be dominated by star formation. The UV and ra- 
dio morphology suggest a star formation rate den- 
sity of 15 — 50M ( r l yr~ 1 kpc~ 2 , within the range of 
locally observed starbursts JSolomon fc Vanden Boufj 
l2005tlSanders fc Mirabe]|ll996l ). However, without a high 
resolution map of the mm emission and gas it is not pos- 
sible to form a clear picture of how the star formation is 
distributed. No X-ray flux is detected in a 200ks Chan- 
dra exposure placing the X-ray to FIR luminosity ratio 
in the star f ormation dominated regime (| Alexander et all 
12005112003ft . and the ra dio to FIR flux r atio falls on the 
local starburst relation (lYun et alJ l2001). However, the 
limit on the X-ray to radio luminosity ratio does not rule 
out an AGN . The SFR inferred by the diffuse UV and 
optical emission can explain the FIR emission if the two 
are spatially related. Finally, the Ly-a line is narrow and 



AGN emission lines such as broad C IV or N V are not 
observed in the optical spectra, so any AGN must either 
be heavily obscured and/or outside the spectrograph slit. 

4. IMPLICATIONS FOR GALAXY FORMATION 

This object is a likely progenitor for the massive (> 
10 11 M Q ) old (> 2Gy r) early type systems s een in large 
numbers at z ~ 2 dMcGrath et al.l [20071: iKong et al.l 
120061: iDaddi et al.1 120051: ICimatti et al.ll200l . The mor- 
phology and spectral properties of the passive galax- 
ies indicate they fo r med in a single bur s t at z > 4 
(ICimatti et al.l 120081 : iStockton etafl 120081 : IDaddi et al.1 
12005ft . However, the density of passive z ~ 2 system s 
is ~ 10- 4 Mj?c- 3 (IKong et al.ll200l IDaddi et al.ll2005ft . 
which is too high to be ex plained by the previ ous mm 
source redshift distribution (|Cimatti et al.ll2008ft . 

The discovery of this object and other recent studies 
suggest the fraction of z > 4 mm sources may be higher 
than previously thought. An 850^m flux limited sample 
is equivalent to a star formation rate limited sample at 
0.5 < z < 7.5 and the z ~ 2 objects must have formed 
by z ~ 4 in order to have sufficient time to evolve into 
passive systems. With this redshift range and a star for- 
mation duration (duty cycle) of 50Myr, a surface density 
of ~ 200 objects per square degree is required for these 
sources to be progenitors of the z ~ 2 passive galaxy 
population. Objects brighter than 4mJy at 850/iin would 
have a sufficiently high SFR to form > 10 11 M Q passive 
systems within 50Myr and the density of such sources is 
sufficient to form the z ~ 2 passive gal axies if ~ 30% 
of them are at z > 4 (jBorvs et al.l [2003). This fraction 
is well within the range recent studies place at z > 4 
(I Younger et all 120071 : IWang et alJ l2006t IChapman et~aTl 
2005^ 
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